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ABSTRACT 

We study by numerical methods a particular kind of SU(N) Yang-Mills solutions of the 
Euclidean equations of motion which appear on the torus when twisted boundary condi- 
tions are imposed. These are instanton-like configurations with the peculiarity of having 
fractional topological charge. We focus on those solutions with minimal non-trivial action 
S = Stt^/N and extract their properties in a few different cases, paying special attention 
to the iV — s> oo limit. 
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1 Introduction 



Quite some time ago 't Hooft IQ pointed out some very special features which arise while 
formulating gauge theories on a torus (a review can be found in reference 0). They are 
related to the freedom in choosing boundary conditions for the gauge potential: since 
only local, gauge invariant quantities are required to be periodic, periodicity of the gauge 
potential has to be satisfied only up to a gauge transformation. Under displacement of the 
gauge potential by a torus period li, in direction z> 

A^(x + /,z>) = f2,(a;) A^(x)f2t (x) - (3.) , ^ = 1, . . . , 4, (1) 

with fl^{x^^i^) SU(N) matrices, also known as the twist matrices. The choice of such 
matrices is arbitrary up to a consistency condition derived from the requirement of single- 
valuedness of the gauge potential. Built in two ways from A^(xi^, Xp) (through Afj_{xy-\4^, Xp) 
and A^{xy,Xp + lp)) single-valuedness for Ap{xy + ly,Xp + lp) implies 

np{x, + iu)n,{xp) = n,{xp + ip)np{x,) Exp (^_!^![^^ (2) 

with Tipjj a gauge invariant antisymmetric tensor of integers, defined modulo N and in- 
dependent of X (this twist factor is allowed due to the invariance of under a gauge 
transformation with an element of the center of SU(N)). Indeed, the actual choice of 
the twist matrices is irrelevant and only the consistency conditions given by n^j/ matter. 
Whenever 7^ (mod N), for some /i, we say the boundary conditions are twisted. 
The twist is reflected in a gauge invariant way through the non-trivial periodicity of the 
Polyakov loops, defined on the torus as 

Cp{x) = ^Tr (L^) = ^Tr (^Texp A.dx''^ f2^(a;') Texp A.dx"^ (3) 

with 7^(a,6) a straight line in the positive direction starting at a and ending at h and x' 
the border of the torus patch. Periodicity holds only up to the twist factors, i.e. 

/:,(x + u) = Exp(-^^)/:,(x). (4) 

With twisted boundary conditions the topological charge Q is no longer necessarily an 
integer: 

Q = J Tr (f^^F^,,) d4X = u - ^, with z/, k G ^, (5) 
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K is associated to the matrix of integers iji^y) through 

1 _ 

= -^n^un^u = k - rh, (6) 

with ki = riQi, riij = eijkrrik. The mathematical proof of this relation for the topological 
charge can be found in 0. Then Q is fractional and proportional to 1/N whenever k-m ^ 
modulo N (non-orthogonal twist). This implies, through Schwarz-inequality, that the 
action of any configuration is bounded from below by 



K 

V 

N 



(7) 



with the bound saturated for self or anti-self dual configurations (F^^, = ±F^,^). It is clear 
that whenever k 7^ (modulo N) there is an obstruction for zero-action configurations. 
Minimal action is attained in such cases if \v — k/N| = 1/N with S = Stt^/N. These are in 
fact the kind of solutions we will describe in this paper. 

Some of these fractional charge solutions have already been found either analytically 
or numerically, 't Hooft has explicitly constructed non-abelian solutions with constant 
field strength which turn out to be (anti-) self-dual whenever the sides of the torus satisfy 
certain relations (see for details or the appendix at the end of this paper). There are 
also a few numerical studies of solutions with non constant field strength. The first one 
is presented in reference and it is obtained there the fractional charge solution with 
\Q\ = 1/2 and S = 47r^ for the SU{2) group, on a L'^ x T torus with T ^ L and satisfying 
twisted boundary conditions given by the twist vectors m = (1, 1, 1) and k = (1, 1, 1). A 
full parametrization of the field strength F^j^, and of the gauge field for this solution 
is presented in reference 0. Another SU(2) solution is presented in |^, in this case the 
fractional charge solution with \Q\ = 1/2 and S = Att"^, on a torus X with T > L 
and satisfying the appropriate boundary conditions to have the properties of a vortex. The 
same kind of solution for the SU(3) group, with \Q\ = 1/3 and S = 87r^/3, is presented in 
reference 0, and the generalization to SU(N) group with > 3 can be found in 0. In 
this article we present a numerical study of SU(N) solutions, with charge \Q\ = 1/N and 
action S = 8n'^/N, and living on a x T torus with T ^ L. Some preliminary results 
have been presented in [ITU . 
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These configurations are interesting by itself from a mathematical point of view, and 
physically interesting for their possible relevance in low energy phenomena like the confine- 
ment property or the breaking of the chiral symmetry. As has been pointed out in ||^ these 
solutions may play a role in the properties of the theory in the limit of large number of 
colors, N. One of the arguments to question the contribution of instantons to long-distance 
phenomena as confinement is based on the large N expansion [|ll[] . Any instanton mediated 
interaction is suppressed by the semi-classical factor exp(— 87r^/(7^), since the large N limit 
is achieved while keeping g'^N fixed, integer charge instantons are (at least in the dilute 
gas picture) naively suppressed by exp(-N). The argument no longer holds for twisted in- 
stantons with action Svr^/N. Another interesting point is the possible relation between the 
center vortex picture of confinement, proposed in |T^ and now being investigated p!3HT7[, 



and fractional charge solutions. As have been pointed out in [0-^], it is possible to built 
vortex configurations in from solutions of the Yang Mills equations of motion in T^. We 
also want to mention the model of confinement based in fractional charge objects presented 
in reference [18|, and some favourable results shown in |T9[| . 

The paper is structured as follows. In section two the numerical method to obtain the 
solutions is briefly described. We will be interested in solutions living on a x T torus 
with T ^ L which, in the limit T —>■ oo, represent vacuum to vacuum tunneling. The 
analysis will be restricted to spatial twist m = (1, 1, 1). Section three presents a detailed 
analysis of these solutions. Our conclusions are presented in section four. Finally, we 
include an appendix with the analytic solutions obtained through 't Hooft construction. 
Their relation with the cases we have studied is discussed through the text. 



2 Numerical minimization of the action 

To generate numerically the minimal action configurations we follow the method which 
has allowed to successfully extract these kind of solutions for other sizes of the torus and 
values of the number of colors in references P-[TO[|. We use the standard discretization 



of Yang-Mills theories on the lattice |2^. We work on A^^ x Nt lattices, Nt 3> Ng, with 
variables defined on each link of the lattice taking values on N x N unitary matrices U^(n). 
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The lattice action used is the Wilson action, 

5^ = Tr (l - U^(n)U,(n + /i)Ut (n + i))ijl{n)) , (8) 

n,fi,u 

where /i and u specify directions (from 1 to 4) and /t , are unit vectors along the 
corresponding direction, = 1, ...,N^. 

The link variables U^(^) satisfy the (twisted) boundary conditions, 

V^{n + N,u) = n,{n)lJ^{n)nl{n + il)) , (9) 

where N4 = Nt, Ni = Ng, i = 1,2,3 and are the twist matrices with consistency 
condition, 

n^{n + N^u) n,{n) = n,{n + N^fi) Vt^{n) Exp(-27rzn^,/N) . (10) 
It is possible to make a change of variables 

U^(Ar^, n,) = U^(Ar^, n,)fi^(n,) U^(n^ ^ A^^, n,) = U^(n^ ^ N^, n,) (11) 
such that the new link variables are strictly periodic. In terms of the new links 

5^ = 5: Tr (1 - Z;^in) V^in) U.(n + /i) ^^71 + 0) \Jl{n)) , (12) 

where Z^j^yin) G take the values: Z^u{n) = 1 for all plaquettes except the one at the 
top-right corner in the {fi, v) plane which is equal to Exp(— 27rm^i,/N). 

The strategy to obtain the solution is minimize the lattice action with respect to the 
variable \J^{n) (this minimization procedure is usually known as cooling). We use the 
Cabibbo-Marinari-Okawa algorithm |21|] in which each link variable is updated in the way: 
U^(^) ~^ AU^(n), where A is a SU(N) matrix built from a SU(2) matrix a which is 
embedded into one of the A^(A^ — l)/2 subgroups of SU(N). Once we obtain the matrix a 
minimizing the new action S'iy(AU^(n)), we update the link variable U^(n), and repeat 
the procedure for all the N{N — l)/2 subgroups of SU(N) and for all lattice sites. This 
constitutes one cooling sweep. We iterate this procedure up to we obtain that the Wilson 
action is stable with a given precision (in this work, the eight relevant digit) and close to 
the value of the expected continuum action: S = 8tt'^/N. 
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3 The solutions 



As mentioned in the introduction, we are interested in solutions with minimal non-trivial 
action, 



N ' 



(13) 



on a volume [-L/2,L/2f x [-T/2,T/2], with T :$> L. When T ^ oo these solutions 
represent vacuum to vacuum tunneling. 

We have restricted our analysis to the following non-orthogonal twist tensors: 

1. Spatial twist, always fh = (1, 1, 1). 

2. Temporal twist, two cases: 

— * 

• k — (1,0,0) for N — 3,4,5,8,10 the solution is in this case anti-self-dual, 
Q = -l/N. 

• k = {n, n, n) for N = 3n + 1 = 4, 7, 10, 13, 19, 25 the solution is here self-dual, 
Q = l/N. 

A list of all the lattices analyzed is presented in Table 1. 

From the lattice configurations wc can easily derive information concerning continuum 
quantities. Part of it can be extracted in a gauge invariant way, such is the case for 
instance of the eigenvalues of the field strength or the Polyakov loops. However to derive 
information concerning the gauge potential gauge fixing is needed. 

The continuum field strength tensor is extracted, up to O(a^), from the clover average 
of the plaquette: 

Qnu{n) = - 



(14) 



through 



1 1 
a^2i 



Q,An) ^ %M - -Tr (Q^,(n) - Qt,(n)) 



In terms of the gauge fixed links the gauge potential is. 



[n+-)a 



1 1 

a 2i 



U^/(^) - K^Hn) - ^Tr (U^^(n) - V^^^n 



N' 



(15) 



(16) 
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The Polyakov loops L^(x) are simply given on the lattice by the ordered product of the 
/i-links corresponding to the path 7;x(x) in Eq. (||). 

Non-gauge invariant information about the configurations will be presented in the tem- 
poral gauge: A4 = 0. In addition we fix Aj(t = — cxd) = which is allowed because in the 
T — s> 00 limit fractional instantons describe vacuum to vacuum tunneling. In this gauge 
Aj(t = 00) = —tfl4difl\, with the temporal twist matrix, and the spatial twist matrices 
are constant. This is not yet a complete gauge fixing, we still have the freedom to make a 
global gauge transformation and also to multiply the twist matrices by an element of the 
center of the group. We have made use of the global gauge transformation to bring the 
spatial twist matrices to a particular form. An explicit construction of constant spatial 
twist matrices compatible with the twist rh = (1, 1, 1) can be easily found, following 't 
Hooft i: 

rjg = Q = P^"^ ^^1 = PQ^-^ 

with p an integer number taking values p = 1,2, ...N, and P, Q the matrices. 



(17) 



/ 











where 



... 1 

V (-1)^+1 ... 0/ 

exp{i2Tm/N) with n = 0,1 



Q 



J.tt{1-N)/N 



/0O 





V 



\ 




;i8) 



N- 1. 

The invariance under multiplication by an element of the center of the group is fixed 
by imposing that the Polyakov loops take the value Ae^"^ at the position where the energy 
density of the solution is maximal. 

On the lattice the A4 = 0, Aj(t = — T/2) = gauge is implemented by transforming 
the corresponding link variables to the identity. The gauge transformation, uj{n), which 
implements the change, is constructed in the following way: choose a point in the time slice 
if: = —T/2, i.e. n° = {rit = u!{n) is the product of the link variables along a certain 

path connecting with n = {nt,nx,ny,nz). In particular we choose n° = (1, 1, 1, 1) and 
the path such that it reaches the point n first in the x direction up to n^, then in the y 
direction up to Uy, in the z direction up to and finally in the t direction up to n^. In 
this gauge the information about the twist matrices is encoded in the links Uo(^t = Nt), 



Table 1: Set of studied solutions. All quantities are defined on equation In the first column we label 
the solutions for reference in other tables. 



Sol. 


Group 


— * 

k 


Size 


5N 
87r2 


Q N 


87r2 


•S'tN 
87r2 




SU(3) 


(1,0,0) 


6^ X 18 


0.95409 


-0.95370 


0.47817 


0.47592 




SU(3) 


(1,0,0) 


73 X 21 


0.96623 


-0.96602 


0.48393 


0.48230 




SU(3) 


(1,0,0) 


8^ X 24 


0.97413 


-0.97401 


0.48768 


0.48645 


I.l 


SU(3) 


(1,0,0) 


10^ X 30 


0.98343 


-0.98338 


0.49210 


0.49133 


1.2 


SU(4) 


(1,0,0) 


8^ X 32 


0.98273 


-0.98267 


0.49168 


0.49105 


1.3 


SU(5) 


(1,0,0) 


8^ X 40 


0.98742 


-0.98738 


0.49387 


0.49355 


1.4 


SU(8) 


(1,0,0) 


53 X 40 


0.98377 


-0.98368 


0.49194 


0.49183 


1.5 


SU(IO) 


(1,0,0) 


6^ X 48 


0.99204 


-0.99201 


0.49465 


0.49739 


II.l 


SU(4) 


(1,1,1) 


8^ X 32 


0.98275 


0.98268 


0.49166 


0.49109 




SU(7) 


(2,2,2) 


33 X 21 


0.94440 


0.94346 


0.47284 


0.47155 




SU(7) 


(2,2,2) 


53 X 35 


0.98020 


0.98008 


0.49026 


0.48994 


II.2 


SU(7) 


(2,2,2) 


8^ X 56 


0.99229 


0.99228 


0.49613 


0.49617 


II.3 


SU(IO) 


(3,3,3) 


4^ X 40 


0.98192 


0.98180 


0.49100 


0.49092 


II.4 


SU(13) 


(4,4,4) 


43 X 52 


0.98801 


0.98794 


0.49399 


0.49401 




SU(19) 


(6,6,6) 


2^ X 38 


0.97324 


0.97200 


0.48729 


0.48594 


II.5 


SU(19) 


(6,6,6) 


3=^ X 57 


0.98812 


0.98804 


0.49405 


0.49406 


II.6 


SU(25) 


(8,8,8) 


2^ X 50 


0.98205 


0.98187 


0.49109 


0.49096 



Ui(nj = l,nj = A'^s), the latter are rotated to the form indicated in Eq. ([T7| ) . 
3.1 Gauge-invariant quantities 

1. Global quantities. In Table 1 we give the values obtained for the action 5, electric 
and magnetic parts of the action, and respectively, and topological charge Q: 

S = \Jtt (F^,F^,) d,x = Jtt + B^) d,x 

Se = Jtt (e^) d^x Sb = Jtt (b^) d^x 

1 /" ~ 1 /" 

Q = J Tr (f^,F^,) d^x=—J Tr (eb) d^x (19) 

where Ej = F4j and ^eijkFij = B^. We can see from the data that the configurations 
obtained are (anti-) self-dual to a very good degree, being therefore solutions of the Eu- 
clidean equations of motion. Those values are very near to the continuum expected values 
^N/87r2 = 1, 5eN/87r2 = 0.5, Sb^/Sir^ = 0.5 and QN = ±1. 
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The first thing we should check is the scahng of the solutions. Without loss of generahty 
we set the spatial length to Is — 1, being then the lattice spacing a — l/Ns- To see how 
the continuum limit a — >^ is approached we vary the lattice spacing while keeping all 
other parameters fixed (among them the ratio Nt/Ns). We fit the N and dependence 
of the action to the expression 5'N/87r^ = 1 — Aa^/(N\/N) and obtain that for the value 
A = 8.893 the data in table 1 are well described with errors smaller than the 0.2% . From 
this fit we understand how the continuum limit is approached for any value of N, and 
also that the N dependence is such that the lattice corrections decrease with increasing N. 
This property will be discussed further later on, it implies that for large N we can obtain 
good continuum results already on rather coarse lattices, this is a rather general property 
which, as we will see, affects other quantities apart from the integrated action and charge 
densities. 

2. Energy profile. The energy profile, defined as 

e{t) = y Tr (Ef{x, t) + B^{x, t)) d^x, (20) 

is located on a region of size ~ and has only one maximum for all values of N up 
to = 13 (instanton profile) and a double peak structure for the values = 19, 25. In 
Figures la and lb we show the scaling with the lattice spacing for the sohitions with = 3 
and N — 19. We can see that points coming from lattices with different sizes describe very 
similar curves, scaling towards the same continuum function. 
For values of N up to A" = 13 , e(t) is well fitted by 

'^^^^ ^ A cosh{wt) + Bt^ + C ■ ^^^^ 
The values obtained for the parameters A,B,C,w are given in table 2. For the values 
N = 19, 25 we fit to the expression (0(t — to) + 4>{t + ^o))/2 and also in table 2 we give the 
results of these fits. 

The N dependence is such that 

^ ^ (22, 

as illustrated in figures Ic and Id where N^e(t) is plotted as a function of t/N for k = (1, 0, 0) 
and N = 4, 5, 8, 10 in figure Ic and for k — (n, n, n) and N = 4, 7, 10, 13, 19, 25 in figure 
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Figure a 



Figure b 



N=3 



Squares --> Ns= 6 
Circles --> Ns= 7 
Diamonds --> Ns= 8 
Ellipses --> Ns=10 




I 

262.5 
175 



Squares --> Ns= 
Circles --> Ns= 2 
Ellipses --> Ns= 3 




Figure d 





A\ Circles --> N=7 




Squares --> N=10 




/ \ Diamonds --> N=13 




f'^*\ Empty ellipses --> N=19 




/ i Ellipses --> N=25 




* 1 




^ 


/ \ 


\ 

Y^'^r.v, 



t/N 



Figure 1: In figures a and b it is shown the energy profile as a function of time calculated with different 
lattice sizes for the solutions with N=3 and N=19, respectively. In figures c and d it is shown the N 
behaviour of the energy profile, plotting iV^ times the value of the profile as a function of t/N , in figure c 
for solutions with k = (1, 0, 0) and in figure d with k = (n, n, n). The continuum curves are the fits to the 
function Ell given in table 2. 
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Table 2: Results of the fit to the energy profile. 



Sol. 


N 


A 


B 


C 


N w 


to 




I.l 


3 


0.00863 


0.03009 


0.02132 


17.95 


0.0 


0.0001500 


1.2 


4 


0.02371 


-0.01394 


0.02588 


17.10 


0.0 


0.0000468 


1.3 


5 


0.04700 


-0.04251 


0.02766 


16.51 


0.0 


0.0000490 


1.4 


8 


0.11399 


-0.07989 


0.07554 


17.54 


0.0 


0.0000016 


1.5 


10 


0.22671 


-0.20144 


0.08006 


17.26 


0.0 


0.0000074 


II.l 


4 


0.02357 


-0.01222 


0.02602 


17.11 


0.0 


0.0000454 


II.2 


7 


0.07820 


-0.03956 


0.06456 


17.63 


0.0 


0.0000037 


II.3 


10 


0.17808 


-0.15184 


0.12957 


18.23 


0.0 


0.0000131 


II.4 


13 


0.41834 


-0.34723 


0.13673 


17.94 


0.0 


0.0000427 


II.5 


19 


1.79852 


-0.59304 


-1.08299 


18.93 


1.24 


0.0000542 


II.6 


25 


2.02385 


-0.59827 


-0.82069 


22.86 


1.74 


0.0000168 



Id. We also plot in figure Id the energy profile for the abelian solution described in the 
appendix (equations and We plot N'^t{t) = 247r^ for the values oit/N between 

t/N = —1/6 and t/N = 1/6. This is the profile for the self dual abelian solution in the 
N — > oo limit. We can see that the energy profiles of the solutions with k = {n,n,n) are 
approaching the one of the abelian solution in the N — > oo limit. 

3. Action density. Defined as, 

Six,t) = TT{E^{x,t) + B^ix,t)) . (23) 

For values of N up to = 13, the action density has only one maximum what we will 
call the center of the instanton. We fit the center and their first nearest neighbours to the 
expression, 

Six, t) = So(l-Y. ^^^^^ - ^-^^] (24) 

where 5*0 is the height, x^,t^ the position and ai,at the width of the maximum. The values 
obtained are shown in table 3. For the values A^ = 19, 25 we observe two maximum in the 
action density and we make the same fit for each one. The results are also shown in table 
3. In both cases, the errors are obtained from the difference between the same quantities 
calculated for the electric and magnetic part of the action. We can see that with increasing 
N all maximum become spatially flat. In fact, this is a general and important property of 
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Table 3: Sizes and heights of the maximum in the action density. 



Sol. 


N 




at 


So 


I.l 


3 


0.531(4) 


0.389(2) 


63.0(2) 


1.2 


1 


0.760(6) 


0.181(5) 


27.56(12) 


1.3 


5 


1.03(3) 


0.583(4) 


15.68(12) 


1.4 


8 


2.90(4) 


0.980(4) 


5.38(2) 


1.5 


10 


6.0(2) 


1.535(5) 


3.25(1) 


11.1 


4 


0.763(10) 


0.480(5) 


27.65(12) 


II.2 


7 


1.98(3) 


0.814(3) 


7.29(4) 


II.3 


10 


6.53(7) 


1.51(1) 


3.253(7) 


11.4 


13 


19.5(1.5) 


5.2(8) 


1.7835(7) 


11.5 


19 


290(20) 


2.3(2) 


0.794(8) 


11.6 


25 




2.4(2) 


0.448(3) 



the solutions; when N is large some quantities, among them the action density, are spatially 
independent. In particular this implies that all the coordinate dependence of the action 
density comes through the time dependence of the energy profile e{t) defined above. This 
fact allows to easily understand the decrease of the lattice artifacts with increasing N since 
generally constant fields give rise to a much smoother, continuum-like, behaviour. 

4. Eigenvalues of F^j,. 

Since the solution is (anti) self-dual we only give the results for Bj. The main properties 
for the eigenvalues of F^^ are illustrated in figures 2a, 2b, 2c and 2d. We only show the 

— * 

results for the solutions with twist k — {n, n, n) because the same properties are obtained 

— * 

for the solutions with twist /c = (1, 0, 0). 

— * 

In figure 2a we show the eigenvalues of Bi for the solutions with N — 7 and A; = (2, 2, 2). 
Very similar results are obtained if we plot B2 or B3 instead of Bi. That we plot is the 
spatial average of each eigenvalue as a function of time. The error bars mean spatial 
dispersion of the eigenvalues (difference between the maximum and minimum value of the 
eigenvalue at each temporal point). The first property we observe is spatial independence 
of the eigenvalues. This property also holds for values of > 7 as shown in figures 2b, 2c 
and 2d. Prom figure 2a we observe that we obtain good results from very coarse lattices. 
In this figure we plot points coming from lattices with the following sizes: 3^ x 21, 5^ x 35 
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and 8^ x 56, being the results almost independent of the lattice size. This property also 
holds for bigger values of A^. 

The behavior with N is shown in Figure 2b in which we plot the first, second and 
N*'' eigenvalues of N x Bi as a function of t/N for N = 7, 13,25 and k = {n,n,n). The 
meaning of the points and the error bars is the same as in figure 2a. The first eigenvalue 
is approximately independent of N and the other N-1 become degenerate for increasing N. 
This structure is very similar to the one of the selfdual abelian solution described in the 



appendix (equations 51, ^ and [5^) , the first eigenvalue of N x Bi takes the value N x ^ 



Bo = ^(Bi + Ba + Bg) 



Bl = — (-2B1 + B2 + B3) 



2lT 

Ui iV ^ Ul hiilSilS LlitJ VfcllUt; iV A 

and the other — 1 are equal to the value — x ^^^^^ . To compare we plot on figure 
2b the first eigenvalue for the abelian solution in the N —>■ 00 limit, the value 27r for the 
interval — 1/6 < t/N < 1/6. We see that the first eigenvalue of for the solutions with 
k = {n,n,n) are approaching the one of the abelian solution in the N — ^ 00 limit. 

The symmetry properties of the spatial twist vector m suggest us to consider the fol- 
lowing combinations of Bi fields, 

1 

73 
1 

Bt = (-B2 + B3) (25) 

one parallel to the m vector and the other two perpendicular. Also interesting is that, 
if there is a common component in color space for Bi fields, we will see this component 
appearing in Bq and not in i?^ and Bt- 

In figure 2c we show the first, second and N*'* eigenvalues of A^ x Bq as a function of 
t/N for the same solutions appearing in figure 2b (with the same meaning for points and 
error bars). We observe the same structure seen in Bj and the expected property if there 
is one common component in color space, the dominant eigenvalue is bigger than the one 
for Bj. We can also observe that this eigenvalue of A^ x Bq is approaching to the shape of 
the one for the abelian selfdual solution in the N ^ 00 limit, in this case 

In figure 2d we show the first, second, (A^— 1)*^ and A^*'' eigenvalues of A^ x B/^ as a 
function of t/N for the same solutions of figure 2b (also with the same meaning for points 
and error bars). The same results are obtained if we show B7- instead of B^. The eigenvalue 
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Figure 2: In figure a it is shown the eigenvalues Ai of Bi as a function of time, for the solution with 
N = 7 and different lattice sizes. To compare with the other figures Ai is multiplied by N and t divided 
by N. In figure b we plot the eigenvalues 1, 2 and N of Bi multiplied by N as a function of t/N, for three 
values of TV. In figure c the same as in figure b, in this case for Bq. In figure d it is shown the eigenvalues 
1,2, N — 1 and A'' of Bz, for the same solutions as in figure b. 
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Figure 3: In figure a the eigenvalues A12 of i[Bi,B2] multiplied by N'^/'^ are plotted as a function of 
t/N for three values of N . In figure b the quantity defined in equation cos^iaij) for i = 1, j = 2 is 



plotted as a function of t/N for the solutions with k = (n, n, n). 

Figure a Figure b 
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structure is completely different to the ones shown before. In this case the eigenvalues are 
distributed in pairs, each pair with two opposite values. Another interesting property is 
that at t/N = the eigenvalues go to zero very fast for large N, being Bq the only one 
non trivial in this limit. 

5. Colour orientation of the field strength. This is studied by calculating the 
commutators of Bj's fields. We show in Figure 3a the eigenvalues of A^^/^[Bi,B2] as a 
function of t/N for the solutions of Figure 2b,c,d (again the error bars represent spatial 
dispersion). The same results are obtained if we plot the eigenvalues of A^^/^[B2,B3] or 
A^^/^[Bi,B3] instead of A^'^''^[Bi, B2]. Only two eigenvalues are relevant being the other 
N-2 very close to zero. We can see a nice scaling with A^^/^ of the two relevant eigenvalues 
at points with \t\/N > 0.1, but not at the center of the solution in which the approach to 
zero is faster. For these quantities the spatial independence for large N also holds. 

To study the abelian content of F^,^ we calculate the quantity: 
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Whenever cos^ (ajj) = 1 the solution is abehan. In Figure 3b we plot cos^ (012) as a function 
of t/N for the solutions with N = 4, 7, 10, 13, 19, 25 and k = {n, n, n). As in other figures, 
the points mean the spatial average of the quantity and errors the spatial dispersion. For 
large N the solutions become abelian at the instanton center (t=0). cos^ (cdj) shows again 
X independence for large N. 

For the solutions with k = (1, 0, 0) we obtain the same behaviour near the instanton 
center, but for T — > ±00, cos^ {cnj) goes towards a N depending constant. 

6. Polyakov Loops and structure of vacuum. In the gauge we have chosen the 
relationship between the Polyakov loops and the twist matrices is specially clear. We have 
A4 = and hence the temporal Polyakov loop directly provides the twist matrix fl^lx). At 
t = —oo Aj is also fixed to zero and in consequence 

hi{x,t = -00) = (27) 

Since the spatial twist matrices are constant, compatibility with the spatial boundary 
conditions for m= (1,1,1) (see Eq. (^)), implies that N£^(t = — 00) = Tr(L^(t = — 00)) =0. 
From the boundary conditions in the time direction: 

A(x,t = T/2) =exp(^z/c,^) A(x,t = -T/2) , (28) 

it is clear that also at t = 00 the spatial Polyakov loops are zero. To characterize the 
vacuum states between which the configurations interpolate we need thus an additional 
quantity provided by '^Czyx = '^''^(X'z^y^x)- Using the twist matrices given in Eq. (p!7D the 
values for Czyx in a vacuum are, 

JO-zyx = exp(^27rp/N) (29) 

where p takes the values p = 1,...,N. There are therefore N different vacua labeled by 
the value of Czyx- Our solutions interpolate between two of them as can be seen from the 
boundary condition in the time direction for Czyx- 

Czyx{xj,t = T/2) = exp (zY^k—j Czyx{xj,t = -T/2). (30) 
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Table 4: Results of the fits to equations ^2, 34 and 



Sol. 


AT 


A 




B 


w 










10*x^ 


Nf 






I.l 


3 


0.453 


334.1 


0.657 


8.04 


0.0188 


2.057 


15.77 


20.52 


13.898 


1.2 


4 


0.472 


75.41 


0.764 


8.20 


0.0125 


2.093 


15.79 


16.87 


0.5014 


1.3 


5 


0.485 


12.96 


0.846 


8.22 


0.0077 


2.099 


15.64 


15.70 


0.4166 


1.4 


8 


0.504 


0.4252 


0.969 


8.35 


0.0221 


2.095 


15.93 


14.88 


0.0411 


1.5 


10 


0.508 


0.0094 


0.998 


7.98 


0.0214 


2.092 


16.02 


14.60 


0.0375 


II.l 


4 


0.288 


0.2311 


0.763 


8.11 


0.0534 


2.094 


15.82 


16.88 


0.4846 


II.2 


7 


0.150 


0.0017 


0.929 


8.11 


0.0025 


2.099 


15.94 


15.07 


0.0868 


II.3 


10 


0.102 


< 10-6 


0.999 


8.19 


0.0334 


2.093 


16.19 


14.70 


0.0353 


II.4 


13 


0.078 


< 10-6 


1.029 


8.05 


0.0925 


2.095 


16.53 


14.52 


0.0428 


II.5 


19 


0.053 


< 10-6 


1.057 


7.96 


0.2738 


2.103 


17.07 


14.32 


0.0425 


II.6 


25 


0.040 


< 10-6 


1.066 


7.88 


0.4323 


2.123 


17.45 


14.21 


0.0389 



We can parametrize the data obtained for Ci, Cq, Czyx as, 



= Mt) e^^'W exp |^(m x f),| (31) 
Coix) = A e'°° expl^^kr] (32) 



N 

^zyx{x,t) = fzyx{t) C ^y^O (33) 

note that with the gauge fixing condition for the Polyakov loops those functions take the 
values: ao = 7r, ai{t = 0) = vr (r = is the maximum of the solution). 

For Ci we make a fit at each temporal point to the spatial dependence of equation |31|. 



The values of yx^/^s obtained for the solutions with temporal twist vector k = {n,n,n) 
are always smaller than 1.24°, 0.27°, 0.07° and 0.004° for the values = 4, 7, 10 and 13 
respectively. For the solutions with temporal twist vector k = (1,0,0) these values are 
2.7°, 1.28°, 0.81°, 0.13° and 0.05° for = 3,4,5,8 and 10 respectively. We extract the 
values of Ci at the spatial maximum of the solution and make a fit to the expression: 

m = ^cosh{wt/N) (34) 

the values of B and w are given in table 4. 

The functions f^y^, aiy^'"-* and a^^yx'^^ parametrizing C^yx in equation |33| are well fitted 
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by: 

fzyxif) — 1 ~ 



zyx 



N cosh{wzyx t/N) 

TT , . t 



TT , . t 



<^zyxit) = n[^ + tagh{v,y^ — ; 



t^gKvzyx (35) 

the values obtained for A^yx, w^yx, Vzyx ^-re given in table 4. And finally Cq only needs the 
constant A to be fitted, their values are given in table 4. 

We compare our results with the Polyakov loops for the abelian selfdual solution de- 
scribed in the appendix, 

1 ,^ ^ /.27rAr-l i\ ^ /.27r,^ 



where — T/2 < i < 7/2 and —0.5 < Xi < 0.5. We can see that the Polyakov loops in the 

— * 

N ^ oo limit for the solutions with k — {n, n, n) are the same as the ones for the abelian 
selfdual solution. 



3.2 Gauge-dependent quantities 

1. Eigenvalues of Aj. After gauge fixing to the gauge described previously, we calculate 

the eigenvalues of Aj. The main properties of these quantities arc shown in figures 4a, 4b, 
4c and 4d. We only show the results for the solutions with twist k = (n, n, n), because the 
same figures are obtained for the solutions with twist k = (1, 0, 0) changing the sign of the 
gauge field, Aj ^ — Aj. 

In figure 4a we show how the eigenvalues of Ai scale towards the continuum limit. We 

— * 

plot these quantities for the solutions with N — 7, temporal twist vector k — (2, 2, 2) and 
lattice sizes 3^ x 21, 5^ x 35 and 8^ x 56. As in previous figures, points values are the 
spatial average of the eigenvalues and error bars mean spatial dispersion. We can see that 
the discretization errors still are important for the solution with size Ng = 3 but are very 
small for bigger sizes. For these quantities we consider in the following lattice sizes with 
A^. >4. 
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Figure 4: In figure a we show the eigenvalues fxi of Ai as a function of time for the solution with N = 7 
and different lattice sizes. In figure b it is plotted the eigenvalues 1, 2 and A'' of Ai as a function of t/N, 
for three different values of N. In figure c the eigenvalues 1,2 and A'' of Aq are plotted as a function of 
t/N. In figure d the eigenvalues 1, 2, A'' — 1 and N of Al multiplied by \/N are plotted as a function of 
t/N for the same solutions shown in figure b 
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The behavior with N is shown in Figure 4b in which we plot the first, second and N*'* 

— * 

eigenvalues of Ai as a function of t/N for N = 7, 10, 13 and k = (n, n, n) (points and error 
bars have the same meaning as before). Very similar results are obtained if we plot A2 or 
A3 instead of Ai. One of the eigenvalues is approximately independent of N and the other 
N-1 become degenerate and approach zero for increasing N. We also show the biggest 
eigenvalue for the abelian selfdual solution described in section 3 of the appendix in the 
N ^ 00 limit. In this limit the function describing this eigenvalue is 27r(-^ + |) for the 
values -I < jf < l- We can see that the first eigenvalue of Ai is approaching the one of 
the abelian selfdual solution in the N ^ 00 limit. 

As for Bj fields, we consider the following combinations of Aj fields, 

Ao = (Ai + A2 + A3) 
Al = (-2Ai + A2 + A3) 

At = ^ (-A2 + A3) . (36) 

In figure 4c we show the first, second and A^*^ eigenvalues of Aq for the same solutions 
appearing in figure 4b (points and error bars have the same meaning as before). We observe 
that the eigenvalue structure is the same one seen for Ai and the expected result if there 
is a common component in color space for Ai, A2 and A3, the first eigenvalue of Aq is 
bigger than the first one for Aj. We also show the first eigenvalue of Aq for the abehan 
selfdual solution in the N — > 00 limit, in this case the function 27r-\/3(^ + |) for points 
— I < ;^ < |- The first eigenvalue of Aq is approaching the one of the abelian selfdual 
sohition in the N ^ 00 limit. 

In figure 4d we show the first, second, {N — 1)*^ and iV*^ eigenvalues of \/N x A^, for 
the same solutions appearing in figure 4b (points and error bars have the same meaning 
as before). Very similar results are obtained if we plot the eigenvalues of A^- instead of 
the ones for A^,. The eigenvalue structure is completely different to the one shown for Aq, 
Aj. The eigenvalues are distributed in pairs, each pair with two opposite values. As can 
be seen from the figure these eigenvalues goes to zero as 1/y/N. This means that A^ and 
At go to zero for large N while Aq is independent of N. 
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Figure 5: In this figure we plot the eigenvalues /X12 of i[Ai, A2] multiplied by iV^/^ as a function of t/N 
for three values of N. 
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2. Colour orientation of the potential. As for the field strength, this is studied 
by calculating the commutators of Aj's fields. We show in Figure 5 the eigenvalues of 
A'"^/^[Ai, A2] as a function of t/N for the solutions of Figures 4b, 4c and 4d (again the error 
bars represent spatial dispersion). Only two eigenvalues are relevant being the other N-2 
very close to zero and it is also clear the scaling with iV^/^ of the two relevant eigenvalues. 
From this figure we conclude that the gauge field Aj become abelian in the large N limit, 
compatible with the properties presented before, Aq is the remaining component in this 
hmit while A^,, A^- goes to zero with N'^/'^. 

3. Twist matrix f24. The temporal twist matrix behaves very differently depending of 
the temporal twist used. Note that their trace is the temporal Polyakov loop given before. 
Wc study their eigenvalues and obtain the following. For temporal twist k = (n, n, n) we 
fit the eigenvalues to the expression: 
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A, = Exp I ^^1^^ + + ^iVTr I j = 2(iV-l)/3 + l,...,iV-l 

An = Exp I — h ^7^ I (37) 

This expression was obtained from the numerical data for the eigenvalues of the twist matrix 
n^. This is a good parametrization for all values of N studied (A^ = 4, 7, 10, 13, 19, 25) 
and the fits to this expression are better for bigger values of N. The interesting point is 
that this expression for the eigenvalues is the one for the eigenvalues of the twist matrix 
O4 given in Q for the abelian solution described in the appendix. 

For temporal twist k = (1, 0, 0) we obtain for large N that the eigenvalues only depend 
on the Xi coordinate, but we have not found a good parametrization in this case. 



4 Conclusions 

In this paper we have presented a set of solutions of the SU(N) Yang Mills equations of 
motion. These solutions are selfdual or antiselfdual, have fractional topological charge 
Q = 1/N and live on the four dimensional torus, T^. We have studied the case when 
the lengths of the torus are x T with T ^ L and with twist vectors m = (1, 1, 1) and 
/c = ^(1, 1, 1) , = (1, 0, 0). Now we summarize the main results we have obtained. 

The obtained results show a clear tendency to describe continuum functions, being the 
different lattice sizes used enough to observe independence of the number of lattice points. 
This property indicates that the obtained configurations describe continuum Yang-Mills 
fields. 

For each value of and twist we always obtain the same solution up to a gauge 
transformation and a spatial translation. This means that we can repeat the procedure to 
obtain another configuration and the differences observed will be a gauge transformation 
and a spatial translation. 

The main characteristic of the solutions are the following 

• The obtained solutions are selfdual or antiselfdual in all the studied cases. We observe 
numericaly that this property is satisfied with a very high precision. This guarantees 
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that these configurations are solutions of the equations of motion. 

The size of the solutions is approximately By size we understand the length of 
the region in the temporal direction in which the core of the solution is included. We 
can see this property in all quantities shown in section 3. For example, we can look 
at the energy profile and check that the most relevant part is located in a region of 
size ~ A^/3. Also in the same region the eigenvalues of the field strength take their 
maximum values and go to zero out of this zone. We can say similar assertions for 
all quantities calculated in this paper. 

The action density has only one maximum for values of A'^ < 13 and a double peak 
structure for bigger values. The spatial dependence of the action density disappears 
with increasing N being only dependent of the temporal coordinate. 

The orientation of the field strength F^i, in color space is very dependent on the 
value of A^. For smaller values of A^ the different components of F^j, are built from 
different components in color space while for bigger values of A^ the same component 
in color space gives the main contribution to the field strength. This property also 
holds for the gauge field A^. This means that, in the A^ — > oo limit, these solutions 
are abelians. 

In the N ^ oo limit one eigenvalue gives the most important contribution to some 
quantities calculated in this article. In this limit, each component of the field strength 
Ff^i, has one eigenvalue which is approximately A" times bigger than the other A^ — 1 
eigenvalues. This property also holds for each component of the gauge field A^. 

Independence of the temporal twist vector k for some of the calculated quantities. 
This is an A^ independent property which holds when the length in the temporal 
direction goes to oo. This property can be seen, for example, in the field strength 

— * — * 

Fut,; if we change the temporal twist vectors used, ki — (1, 0, 0) and k2 — {n, n, n) we 
obtain that Bj'^ = Bj'^ and Ej'^ = — Ej'^. A similar property is held for the gauge 
field, under the change of the twist vectors we obtain the relation Af^ = — Aj'^. 
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We have not succeeded in finding an analytic expression describing the properties of 
the studied solutions. Obviously, the first requirement for an ansatz prepared to find the 
analytical expression is that this ansatz satisfies the previously shown properties. The 
most promising approach seems to be an ansatz based on the similarity of the solution in 
the large N limit with the abelian solution presented in the appendix. Nevertheless, if the 
solutions in this limit coincide with the abelian solution, something singular must happen 
at points \t/N\ = 1/6, being therefore not ease to use this similarity to find the analytic 
expression. To conclude, we hope that all the numerical data presented will be helpful for 
other attempts to find the analytical expression of the solutions presented in this paper. 

A. Analytic solutions. 

't Hooft construction. The and F^,^ fields are built from a diagonal matrix T, in 
the following way 



where a^,y is an antisymmetric tensor and 1^ the length of the torus in the /j, direction. The 
matrix T has the form. 



being k and I integer numbers {k + I = N). To build the twist matrices we use the P 






(39) 



and Q matrices defined in equation |TB|. 
matrices. 



From these matrices we construct another set of 




(40) 



satisfying the properties. 




(41) 



And the ansatz for the twist matrices is 




(42) 
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where Sf^,tfj,,Ufj, and are arbitrary integer numbers. These matrices must satisfy the 
consistency condition, 

i27m. 



n^{x, + iu)n,{x^) = n,{x^ + i^)n^{x,) Exp (^-^^^j ■ (43) 

This condition imposes the following equations for Sfj,,t^,u^ and 

l{t,s^-Us,)=l^ + ^ + A,^ ; ]{v,u^-v^u,) = -k^ + ^ + B,^ (44) 
where A^^, and 5^,^ are integer numbers. To solve these two equations we give to a^^ and 



n^i, the form, 



„(1) „(2) 

n.. = <j+<) ; «^^ = -f-T' ^^^^ 



and then equations in formula 44 are transformed to. 



which can be solved if the following condition is satisfied, 

<^4'-^ = modk ■ ng)ng)=0 modi (47) 

The n^^J and ra^^J tensors are orthogonal twist tensors for a SU(k) and SU(1) group respec- 
tively. 

The topological charge for these solutions is, 

Q = l^I^kl = - ( fc(i)m(i)- + ^(2)m(2)- - ^(^)m(2) - B^^m^A (48) 
iV 4 N \ k I J 

where we have defined the vectors fcj^"-' = nQi^"'\ rrii^"''^ = eijkUjk^"''' /2, with n = 1,2. 

Some examples. Now we give some examples of solutions built using 't Hooft construc- 
tion, and for the torus lengths used in this article: = ly = h = ^ and It oo. In fact, 
our examples will be given for any value of If. The topological charge in our examples is 



given by equation and the action by the following equation. 



N 2 [It \ k I J 



k (m^'^m^'^^ + m(2)m(2):^ _ 2m^'^m^A \. (49) 
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The minimum value for the action, S — Stt^IQ], is obtained when the solution is selfdual 
or antiselfdual. This condition imposes a value for the temporal length It, obtained solving 
the equation, 

^(1) ^(2) /^(l) ^{2)\ 

the positive sign for selfdual solutions and the negative sign for antiselfdual solutions. Now 
we give some examples for the twist vectors used in this article: 

— * 

1. Solutions for twist vectors m — (1,1,1) and k — {n,n,n), with N — 3n + 1 for 
n — 1,2,3,.... We choose k — N — 1 and I — 1 and the twist vectors in subspaces 
SU(k) and SU(1) as, 

m^') = (1, 1, 1) k^^^ = (n, n, n) ; m^'^ = (0, 0, 0) k^''^ = (0, 0, 0) 

The topological charge for this solution is Q — 1/N and it is selfdual when It — n. 

2. Solutions for twist vectors m = (1, 1, 1) and k = (1, 0, 0). We choose the twist vectors 
in SU {k) and SU {I) subspaces as. 



m 



(1) 



^(^) = (0,0,0) ; m(') = (0,0,0) = (1,0,0) 



This choice works for any values of fc, The topological charge for this solution is 
Q — —1/N and in this case it is not possible to solve the self duality equation. 

Changing the gauge. Now we change the gauge for one of the solutions described before 
to the gauge used in section 3. Wc choose the solution given in example 1. In this case 
A; = - 1, / = 1, nW = n^^ and ng] = 0. The fields are, 

we remember that the torus lengths were = ly = h = and It can take any value, and 
the twist vectors were m = (1, 1, 1) and k = {n,n,n), with iV = 3n + 1. The action and 
the topological charge take the values, 

^~ N 2[n^lt) ^~ N 
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the twist matrices are, 

«„w=p?Q;«i5.p{-^^i:!^T) (52) 

and the values for s^, i^, 

s^ = Q S2 = l S3 = N-2 S4 = 2(iV - l)/3 

The gauge used in section 3 is: A4 = and Aj(t = —00) = 0. For our example we can not 
use the same gauge because t is finite and F^,^ ^ for any value of t (this is the condition 
needed to put Ai{t) = at some point t). The most similar gauge condition is the following 
one, 

A4 = 0, A3(i = 0) = 0, A2(i = 0,z = 0) = 0, Ai(i = 0, z = 0, y = 0) = 

because all links associated with these fields were rotated to the identity. The field A^ and 
the twist matrices in this gauge are. 



A, 



4 



A3 = ^^TTT -t 








27r 




N{N - 


1) 


27r 




A^(A^- 


1) 


27r 




N{N - 


1) 



Ai = tttI^t (^t + ^ - y) (53) 

= P^'^Q^'^ Exp (y - •^■)} 

= P,^iQi*i (54) 

We compare along section 3 the results obtained for large values of N with this solution 
(with It — n) because some properties are very similar for both solutions in this limit. 
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Finally, we give the value of the Polyakov loops for this example, 



note that this quantity is gauge invariant and could be calculated with the two different 
given before, obtaining the same result. 
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